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Abstract 

Triangulation and trilateration data from two geodetic networks located between the 
western edge of the Great Valley and the San Andreas fault have been used to calculate 
shear strain rates in the Diablo Range and to estimate the slip rate along the Calaveras and 
Paicines faults in central California. Within the Diablo Range the average shear strain rate 
was determined for the time period between 1962 and 1982 to be 0.15± 0.08 prad/yr, with 
the orientation of the most compressive strain at N16°E ± 14*. The orientation of the 
principal compressive strain predicted from the azimuth of the major structures in the region 
is N25°E. We infer that the measured strain is due to compression across the folds of this 
area; the average shear straining corresponds to a relative shortening rate of 4.5 ± 2.4 
mm/yr. From an examination of wellbore breakout orientations and the azimuths of P-axes 
from earthquake focal mechanisms the inferred orientation of maximum compressive stress 
was found to be similar to the direction of maximum compressive strain implied by the 
trend of local fold structures. Our results do not support the hypothesis of uniform fault- 
normal compression within the Coast Ranges. From trilateration measurements made 
between 1972 and 1987 on lines that are within 10 km of the San Andreas fault, a slip rate 
of 10-12 mm/yr was calculated for the Calaveras-Paicines fault south of Hollister. The slip 
rate on the Paicines fault decreases to 4 mm/yr near Bitter. 
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Introduction 

Although most of the relative motion between the Pacific and North American plates in 
central California is accommodated by slip along the San Andreas fault, distributed 
compression and right-lateral strike-slip motion on faults subparallel to the San Andreas 
fault also occur between the western edge of the Great Valley and the continental 
escarpment [Gawthrop, 1977; Page, 1981; Crouch et al., 1984; Minster and Jordan, 1984, 
1987; Eaton, 1984; Namson and Davis, 1987], Compression across the Coast Ranges is 
thought to be oriented nearly perpendicular to the trend of the San Andreas fault and has 
been proposed to result from a combination of slightly convergent relative motion between 
the Pacific and North American plates and low shear strength along the San Andreas fault 
[Zoback et al., 1987]. The majority of recent earthquakes in the Coast Ranges occur on the 
San Andreas fault (Figure 1). There is moderate, diffuse activity in the Diablo Range, 
located between the Great Valley and the San Andreas fault. The May 1983 occurrence of 
a large earthquake (M s =6.7) beneath a young fold [Stein and King, 1984; Stein, 1985] 

near Coalinga, California, has focused attention on the importance of understanding the rate 
of deformation to the east of the San Andreas fault in this region and its relation to the 
overall deformation in the Coast Ranges. In this paper we determine rates of deformation 
in the Diablo Range north of Coalinga from a triangulation and trilateration network and 
from line-length changes determined from short-range trilateration measurements within 10 
km of the San Andreas fault. The geodetic data are used to place bounds on the rate of 
deformation across the Diablo Range and to estimate the rate of slip on the Calaveras and 
Paicines faults. These geodetic results are then compared with other geological and 
geophysical data from the Coast Ranges . 
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Tectonic Setting 

In this section we give an overview of the tectonics of the central Coast Ranges. The 
rate and orientation of slip along the San Andreas fault, as well as the magnitude of shear 
stress in the fault zone, have been hypothesized to influence deformation throughout the 
Coast Ranges. We first summarize results from geologic and geodetic studies on the 
segment of the San Andreas fault between Hollister and Coalinga. We then describe the 
geological setting and seismicity of the study region. For comparative purposes, we also 
discuss deformation to the west of the San Andreas fault in central California. 

To characterize horizontal deformation across the San Andreas fault in central 
California, data from near-fault alinement arrays, creepmeters, and trilateration 
measurements at short and intermediate distances have been utilized [Savage andBurford, 
1973; Thatcher, 1979b; Burford and Harsh, 1980; Lisowski and Prescott, 1981]. On the 
northwestern segment (from the latitude of Hollister south to Hepsedam, Figure 2), the rate 
of surface slip (creep) across the San Andreas fault increases from ~ 13 to 32 mm/yr and 
surface slip on the Calaveras-Paicines fault decreases from ~17 to 0 mm/yr. On the central 
section (from Hepsedam southward to the latitude of Coalinga, Figure 2), near-fault and 
intermediate-scale geodetic measurements of right-lateral slip are in good agreement 
suggesting creep at a rate of approximately 32 mm/yr. The rate of slip on the central 
segment of the San Andreas fault estimated from Holocene geological data is 34 ± 3 mm/yr 
at an azimuth of N41 # W ± 2* [Sieh and Jahns , 1984; Minster and Jordan, 1984], 
Southeast of the latitude of Coalinga, shallow slip on the San Andreas fault decreases and 
the width of the zone of deformation increases with the transition to a locked segment of the 
San Andreas fault in the Carrizo Plain. 

The principal structural features of the Diablo Range northeast of the San Andreas 
fault between Hollister and Coalinga are shown in Figure 3, [modified from Dibblee, 
1979]. Within 5-10 km of the San Andreas fault the primary geologic structures are related 
to dextral shear on the San Andreas and Calaveras-Paicines faults. The Diablo Range to the 
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east of this region is a broad antiform which trends approximately N65°W and 
encompasses subsidiary fold structures such as the Vallecitos Syncline [Page, 1985]. 
North of Coalinga the range is pierced by the New Idria diapir of serpentine and Franciscan 
rocks. Along the northeast boundary of the study area is the Ortigalita fault, a high-angle 
fault along or near the contact of Franciscan rocks and the Great Valley sequence to the east 
[ Raymond , 1973]. Trenching of the Ortigalita fault zone shows exposed offsets of late 
Pleistocene and Holocene soils, with possibly as much as 5 km of Quaternary right-slip 
displacement [Hart et al., 1986]. The most recent uplift of the Diablo Range began in 
Pliocene time and most likely accelerated in the Pleistocene [Page, 1981; Page and 
Engebretson , 1984]. 

Eight focal mechanisms from the Diablo Range compiled by Eaton [1985] show a 
mixture of thrust, reverse, and strike-slip faulting. The location of three of the larger 
earthquakes with well-determined focal mechanisms, the October 1982 Idria event 
(Ml= 5.5), the May 1983 Coalinga event (Ml= 6.7), and the August 1985 North Kettleman 
Hills event (Ml= 5.5), are given in Figure 1. The focal mechanism determined for the Idria 
earthquake indicates thrust faulting on a plane oriented N72*E or reverse slip on a plane 
oriented N64°W [Eaton, 1985]. The Coalinga and Kettleman Hills events have similar 
focal mechanisms with slip occurring on fault planes oriented at about N53°W as either 
thrusting on a plane dipping shallowly to the southwest or reverse slip on a plane steeply 
dipping to the northeast [Eaton, 1985; J.P.Eaton, personal communication, 1987]. The 
earthquake focal mechanisms and geological structures in the area suggest two primary 
modes of deformation to the northeast of the San Andreas fault; compression normal to the 
major fold structures of the region and right-lateral strike-slip motion on faults such as the 
Calaveras-Paicines and the Ortigalita. 

In the central Coast Ranges to the west of the San Andreas fault focal mechanism 
studies of the Coast Ranges [Gawthrop, 1977; Eaton, 1984; Dehlinger and Bolt, 1987], 
along with extensive geologic mapping [summarized in Crouch et al., 1984; Slemmons, 
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1987], suggest variable modes of deformation. Between the San Andreas and the 
Rinconada fault (Figure 1) is the seismically quiescent Salinian block. The upper crust is 
composed of high-strength granite which is only weakly folded and sparsely faulted 
[Dehlinger and Bolt, 1987; Slemmons, 1987]. Focal mechanisms from the Rinconada fault 
show a mixture of right-lateral shear faulting on northwest striking planes, oblique slip, 
and reverse faulting [Gawthrop, 1977; Dehlinger and Bolt, 1987]. From field mapping of 
recent offsets, D.B. Slemmons (personal communication, 1987) and E.W. Hart (personal 
communication, 1988) suggests that right-lateral strike-slip faulting is the major component 
of displacement along the Rinconada fault Between the Rinconada and the San Gregorio- 
Hosgri faults (Figure 1) is the Coastal Belt province. The upper crust consists of the very 
heterogeneous Franciscan complex. Focal mechanisms from this region indicate 
dominantly oblique reverse faulting along northwest-trending, northeast dipping planes, 
with P-axes oriented N20’-50°E or S20 < ’-50°W [Dehlinger and Bolt, 1987]. The San 
Gregorio-Hosgri fault system extends from the junction of the San Gregorio and San 
Andreas faults near San Francisco to northwest of the Western Transverse Ranges (Figure 
1). Significant late-Quatemary right-lateral strike-slip offsets have been measured on the 
San Gregorio fault [Clark et al., 1984]. South of Monterey Bay this fault branches into 
several splays with some branches showing primarily right-lateral strike-slip displacement 
and others showing east-up reverse faulting [Hamilton, 1987]. Focal mechanisms from the 
San Gregorio-Hosgri fault indicate right-lateral strike-slip faulting and oblique reverse 
faulting with a right-lateral component on a northeast dipping plane [Gawthrop, 1977; 
Eaton, 1984; Dehlinger and Bolt, 1987]. 
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Geodetic Strain Rates 

We make use of two geodetic networks to estimate current rates of deformation to the 
east of the San Andreas fault in central California. The San Benito triangulation and 
trilateration network spans the Paicines fault zone just east of the San Andreas fault and 
then extends 50 km in a southeast-northwest direction across the Diablo Range (Figure 2). 
For interstation visibility, triangulation and trilateration stations are situated on the highest 
points in a region; thus, many of the stations in the S an Benito network are located near 
anticlinal peaks (Figure 3). The San Benito network is geographically well-suited for 
measurement of rates of deformation between the Great Valley and the Calaveras-Paicines 
fault To examine more localized deformation within the zone 10 km to the east of the San 
Andreas fault, short and intermediate range lines from the U.S. Geological Survey (USGS) 
Coalinga trilateration network were utilized. The geodetic results derived from the 
Coalinga trilateration data provide an update to the slip rates estimated for the Calaveras 
and Paicines faults by Lisowski and Prescott [1981]. To distinguish between engineering 
and tensor shear strain, we denote the former by y and quote the units as [trad, whereas the 
latter is denoted by e and is given in units of [is train [Savage , 1983]. Uncertainities in all 
cases are one standard deviation (a). 

San Benito Network 

A triangulation survey of the San Benito network was conducted in 1962 by the 
National Geodetic Survey. Seventy of the seventy-two directions utilized in this study 
were second-order observations (a priori uncertainty is estimated to be a<j = 0.7”, Federal 
Geodetic Control Committee, 1984) and two were third order (Gd = 1.2"). An angle 
measurement is the difference of two direction measurements. The expected uncertainty in 
a second order angle measurement (a a ) is 2(<Td.; or about 1.0". The order of the 
triangulation denotes the measurement precision, which is determined by survey 
procedures and is reflected in the degree to which internal checks of the data are satisfied 
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[Thatcher, 1979a]. The principal internal check is the triangle closure requirement that the 
angles of each triangle sum to 180° plus the known excess due to the Earth's sphericity. 
Triangle closures indicate that the standard deviation of a single angle is approximately 1.1" 
for these data, very close to the expected value. 

In 1982 a trilateration survey of the San Benito network was performed by the USGS. 
Regional trilateration measurements made by the USGS on line lengths of 10-40 km have a 
precision of approximately 0.2-0.3 ppm [Savage and Prescott, 1973], Angles observed in 
the 1962 survey were compared with angles that would have been measured in 1982. 
Azimuths for each station to station pair for the 1982 trilateration survey were determined 
with a variation of coordinates adjustment [Anderson, 1969]. 

From the changes in angles between the 1962 and 1982 surveys we calculate the 

• • • 
horizontal shear strain rate components yj and Y2- In terms of the strain rate tensor, £{j, 

t i • i * • 

Yl= ei i~ £22 and Y2 = £12 + e 21> where the strain rate tensor is referred to a geographic 
coordinate system in which the 1-axis is directed east and the the 2-axis is directed north. 
The strain component Y2 is equal to the decrease induced by strain in the right angle 
between northward- and eastward-directed lines, whereas Yl is equal to the increase in the 
angle between lines directed northwest and northeast. Thus, Yl measures the rate of right- 

lateral shear across a vertical plane striking N45*W (or the rate of left-lateral shear across a 
vertical plane striking N45’E) and Y2 measures the rate of right-lateral shear across a 

vertical plane striking eastward (or the rate of left-lateral shear across a vertical plane 

striking northward). Angle changes between the 1962 and 1982 surveys were used to 
obtain the two shear strain rate components (Yj and Y 2 ) or > equivalently the maximum 

shear strain rate (y, where y 2 = y^ + y 2 2 ) and the orientation of the vertical plane with 
maximum right-lateral shear (\|/) [Frank, 1966; Prescott, 1976; Sauber et al., 1986]. For 
comparison with the fold structures of the Diablo Range, the orientation of the maximum 
compressive strain ((3) is sometimes given instead of the orientation of maximum right- 
lateral shear (\j/). 
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The angle changes associated with the stations Bitter, Hepsedam, and Panoche were 
significantly larger than angle changes from other portions of the network and were, 
therefore, examined separately. Bitter and Hepsedam are located near the Paicines fault 
zone. Using the 1 1 angle measurements that include either of the stations Bitter or 
Hepsedam y = 0.56 ± 0. 16 prad/yr and = N28‘W ±11*. The shear strain rate across the 
approximately 10 km- wide zone northeast of stations Bitter and Hepsedam, assuming 
right-lateral motion on the Paicines fault, implies a rate of slip of 5 ± 2 mm/yr. Additional 
geodetic data from the Coalinga trilateration network relating to deformation across the San 
Andreas and Calaveras-Paicines faults is discussed in the next section. The relatively large 
angle changes around station Panoche are not so easily explained. The 7 angle changes that 
included Panoche y = 0.76 ± 0.27 prad/yr, y = N58*W ± 9* or (3 - N13*W ± 9*. These 
results are not consistent with either shear strain across the Ortigalita fault zone or 
contraction across the major fold structures of the region. We suspect a survey blunder, 
but the sparse measurements in this part of the network do not provide sufficient 
redundancy to identify such errors. Because the strain accummulation in this part of the 
network was significantly different than the strain estimated from other portions of the 
network, angles which included the station Panoche were omitted from further analyses. 
Results from the northern portion of the network which would be sensitive to slip on the 
Ortigalita fault are unfortunately weakened by this omission. 

The strain rate parameters estimated on the basis of 25 angle changes in the central 
portion of the San Benito network, excluding angles to Bitter, Hepsedam or Panoche, were 
y = 0.15 ± 0.08 |irad/yr, (5 = N16*E ± 14* or y = N29*W ± 14*. The standard deviations 
reflect both misfit and data uncertainties due to measurement imprecision. Although there 
is no significant strain at the 95% confidence limit, the orientation of the principal strain 
directions are consistent with the geological structures of the region. The azimuth of the 
least compressive strain (N74°W ± 14’) is close to the trend of the major fold structures of 


the region (N65*W). The shear is also close to the trend of the major strike-slip faults of 
the region. 

To search for significant strain inhomogeneity, the shear strain parameters were 
estimated from spatial subsets of these data, again excluding the Bitter, Hepsedam, and 
Panoche stations (Table 1). There is a trade-off between improving the precision of the 
strain estimate by using a large number of angles and averaging spatial variations as the size 
of the sampled region is increased. In the first case, the data were broken into two distinct 
groups, one set closer to the Great Valley ("east") and one set closer to the San Andreas 
fault ("west"). If the measured rate of shear strain were due to strain accumulation on the 
adjacent San Andreas fault, the predicted rate of shear strain would be higher in the west 
subnet. Alternatively, if the rate of compressive strain was higher across the folds near the 
Great Valley, the east subnet might show a higher shear strain rate. There is no suggestion 
of a significant rate difference between the strain rates determined from the east and west 
data subgroups. The set of 25 angles were also divided into "north" and "south" sets to 
look for any change which might be associated with along-strike variations on the San 
Andreas and Calaveras-Paicines faults. Although the strain results differ in the subnets, 
particularly in the orientation of maximum compression, the rate and orientation of strain in 
the various subregions are not significantly different than the average value determined 
from the complete set of 25 angles. 

If we assume that the strain perpendicular to the fold structures is homogeneous across 
the region between the western edge of the Great Valley and the western extent of the San 
Benito network (average width of 30 km) and that there is no extension parallel to the fold 
structures, the average strain rate corresponds to 4.5 ± 2.4 mm/yr of shortening. From a 
geological reconstruction of the structures between the Great Valley and the San Andreas 
fault at the approximate latitude of Coalinga , Namson and Davis [1987] have inferred that 
1 1 km of late Cenozoic shortening has occurred perpendicular to the San Andreas fault. 
This quantity can be used to estimate an average rate of shortening perpendicular to the San 
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Andreas fault. If active folding commenced 5 m.y. ago, the average rate of shortening has 
been 2.2 mm/yr. If folding began as recently as 2 m.y ago an average rate of shortening of 
5.5 mm/yr is required. 

A future trilateration or GPS survey could provide a more precise estimate of crustal 
deformation across the Diablo Range and may give the resolution necessary to discern 
spatial variations in the deformation field. Additionally instead of determining only the 
shear strain parameters (yi and 75 ) all the components of the rate of tensor strain (ey) could 
be estimated. Using all the stations in the San Benito network an additional trilateration 
survey in 1987 would have provided ~0.034 |istrain/yr accuracy in the principal strains. 

Coalinga Trilateration Network 

We have employed short and intermediate-length lines from the Coalinga trilateration 
network to update previous estimates [Lisowski and Prescott, 1981] of the slip rates on the 
Calaveras and Paicines faults (Figure 2). Trilateration measurements are made periodically 
by the USGS in central California on a regional scale spanning a 20 -km-wide zone centered 
on the San Andreas fault and including several distinct faults, as well as on smaller (1-2 
km) aperture networks that span a single fault (Pionne, Dry Lake and Tully nets). The 
estimated accuracy in the line lengths determined from the short-range measurements is 
equal to 4 mm [Lisowski and Prescott, 1981]. 

For the region within 10 km to the east of the San Andreas fault we assume that the 
deformation is right-lateral shear across the Calaveras-Paicines fault system. Strike-slip 
type motion across these faults would shorten north-south lines and extend east- west lines, 
whereas contraction across structures such as the Paicines syncline would shorten 
northeast-southwest lines. Line length changes in the Pionne short-range network and the 
lines from Bitter to Hepsedam and Brown to Cross are consistent with the hypothesis of 
right-lateral slip on the Calaveras-Paicines fault system. 
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If the deformation in this region occurs as rigid block motion along the Calaveras- 
Paicines fault system, the observed length changes can be converted into a fault-parallel 
displacement rate from the slope of the least-squares linear fit to the interstation length data 
[Prescott and Lisowski, 1983]. The slip rates determined from lines crossing the fault at a 
low angle are given in Table 2 for the Calaveras-Paicines region just south of Hollister and 
for the Paicines fault zone between Bitter and Coalinga. The slip rate determined from the 
station pairs that cross only the San Andreas fault (Cross/Chalone and Chalone/Bitter 
lines) plus the slip rate from pairs that are completely east of the San Andreas fault 
(Browns/Cross and Bitter/Hepsedam) are shown to be approximately equal to the slip rate 
from station pairs that span both fault zones (Brown/Chalone and Chalone/Hepsedam 
lines). The line length changes across the Calaveras-Paicines fault zone are consistent with 
~10 mm/yr of right-lateral slip. For the Paicines fault further south a lower rate of 3-4 
mm/yr of slip is indicated. For comparison, the late Quaternary slip rate estimated on the 
Paicines fault estimated from an offset terrace riser and an offset hill ~10 km south of 
Hollister is 3.5 - 13 mm/yr [Harms et al., 1987]. 

Some of the lines in the Coalinga network showed a significant change in slope 
beginning in mid- 1978. This change corresponded to an increase in the rate of slip on the 
Calaveras-Paicines system and a decrease in slip rate on the San Andreas fault. A similar 
change after 1979 in the inferred slip rates for the Calaveras and San Andreas faults was 
reported from geodetic measurements from the Hollister trilateration network [G. Gu and J. 
C. Savage, personal communication, 1986; Matsu'ura et al., 1986]. The increase in slip 
on the Calaveras fault approximately coincided with the occurrence of the 1979 Coyote 
Lake earthquake (Ml=5.7) on the northern Calaveras fault . 

The average fault slip indicated by the three short-range networks (Figure 2), again 
assuming simple block motion, is given in Table 3. The rate of slip on the Pionne net, 12 
± 2 mm/yr, is similar to the rate estimated from the line between Cross and Brown which 
crosses the Calaveras-Paicines fault zone (10 ± 1 mm/yr). Between the Dry Lake and Tully 
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networks the rate of creep on the San Andreas fault increases from 27 to 32 mm/yr. South 
of the Tully network the rate of slip determined from near-fault data is the same on the San 
Andreas fault as the rate estimated from the longer trilateration lines. The short and 
intermediate range trilateration measurements document a transfer of slip associated with 
the Calaveras-Paicines fault zone to the San Andreas fault. The narrowing of the region 
that accommodates the ~32 mm/yr slip rate also corresponds to a transition from a complex 
multi-stranded portion of the fault system with locked segments that break in periodic 
moderate to large earthquakes to a geometrically simpler segment that accommodates slip 
through creep. 

Comparison of Geodetic Strain Directions with Principal Stress Directions 
and Orientations of Geological Structures 

In this section we compare the geodetically determined incremental strain directions to 
the orientation of the maximum principal stress (ai) estimated from wellbore breakouts and 
the azimuths of P-axes determined from earthquake focal mechanisms. The orientations of 
principal stress directions determined from breakout orientations, earthquake fault plane 
solutions, and the azimuths of major fold structures in central California have been used to 
infer the state of stress on and near the San Andreas fault and to distinguish between 
models relating the formation of fold structures in the Coast Ranges to motion along the 
San Andreas fault [Zoback et ai, 1987; Mount and Suppe, 1987; Namson and Davis, 
1988]. 

Due to extensive drilling for oil east of the San Andreas fault, a large number of wells 
have been available for measurement of stress-induced wellbore breakout orientations. 
Borehole breakouts are caused by unequal stress concentrations around a borehole wall and 
create elongations of the hole in directions perpendicular to the orientation of the maximum 
horizontal stress [Springer, 1987]. Breakouts are generally oriented northwest-southeast, 
indicating a northeast-southwest orientation for the compressive stress, a direction 
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perpendicular to the major fold axes near the San Joaquin Valley and at a 70*-90* angle to 
the San Andreas fault [Springer, 1987; Zoback et al., 1987; Mount and Suppe, 1987], To 
the east of the San Andreas fault the majority of breakout measurements have been made in 
wells near the western edge of the San Joaquin Valley (Figure 4); comparatively fewer 
wellbore breakout orientations have been obtained closer to the San Andreas fault within 
the area of the geodetic measurements reported in this paper. 

The orientation of maximum principal stress may also be inferred from well- 
determined fault plane mechanisms, notably from the 1982 Idria, the 1983 Coalinga, and 
the 1985 North Kettleman Hills earthquakes (Figure 4). For the Coalinga and the North 
Kettlman Hills events the P-axis orientation is ~N37’E [Eaton, 1985; J.P. Eaton, personal 
communication, 1987], similar to the <7i directions inferred from breakout orientations in 
the same region. The azimuth of the P axis for the Idria event is N12°E [Eaton, 1985]. 

As may be seen in Figure 4, the breakout orientations and the P-axis directions are 
very similar to the orientation of maximum compressive strain implied by the trend of local 
fold structures. Along the western edge of the Great Valley and near Coalinga, the 
direction of ai inferred from the breakout orientations and two earthquake focal 
mechanisms and the orientation of maximum compressive strain inferred from the fold 
trends, is N30-50°E. To the northwest of Coalinga, including the region spanned by the 
San Benito network, the direction of maximum compressive strain estimated geodetically 
and the orientation predicted from the azimuth of the local fold structures and the ai 
direction inferred from one wellbore breakout orientation and the P-axis orientation of the 
Idria earthquake is N12-25*E. Mount and Suppe [1987] suggest that the folds in this area 
have been inactive since the Miocene; the rate of deformation determined for the region 
spanned by the San Benito network suggests, however, that ongoing compressive strain is 
being accommodated across these folds. Within 5 km of the San Andreas fault, the <T i 
direction inferred from one wellbore breakout orientation and the direction of maximum 
compressive strain inferred from azimuths of folds is N30*-60'E. 
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Mount and Suppe [1987] and Zoback et al. [1987] have used contempory in situ 
tectonic stress indicators near the San Andreas fault system in central California to show 
that the direction of maximum horizontal compression is nearly perpendicular to the San 
Andreas fault. The observation of uniform fault-normal compression in a 50-100 km wide 
zone spanning the San Andreas fault has been cited as evidence for low shear stress on the 
San Andreas fault and a decoupling of off-fault normal compression and transcurrent 
deformation [Mount and Suppe, 1987]. The principle stress directions inferred from focal 
mechanisms and breakout orientations, particularly at distances greater than 30 km from the 
San Andreas fault, support the observation of fault normal compression in portions of the 
Coast Ranges. Our results suggest examination in detail of principal stress and strain 
orientations within the Diablo Range do not support the hypothesis of uniform fault-normal 
compression in central California. 

There are several factors which may complicate the interpretation of measurements 
used distinguish between models relating the formation of fold structures in the Coast 
Ranges to motion along the San Andreas fault . Firstly, there may be processes that act on 
different length scales which might explain the change in the ai direction and the trend of 
fold structures at different distances from the San Andreas fault. For example, within a 
narrow zone centered on the San Andreas fault, local geological structures may be due to 
geometrical complexities in the fault trace [Segall and Pollard, 1980] or the rheological 
structure of the fault zone [Horns et al., 1985]. Also, the local stress field orientation may 
differ from the regional stress field due to material heterogeniety. Finally, as will be 
discussed in the next section, there is the additional complexity that deformation in the 
Coast Ranges is related, in the kinematic sense, to the overall deformation across the 
Pacific-North American plate boundary. 
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Relation of Deformation East of the San Andreas Fault to the Overall 
Deformation in the Coast Ranges 

In this section we examine the rates and sense of deformation estimated for different 
subregions of the Coast Ranges, and we discuss the relation of these results to the 
accommodation of Pacific-North American plate motion as constrained by global plate 
models and other regional studies. Finally we briefly suggest future geodetic studies which 
may provide additional constraints on the rate of deformation in the Coast Ranges. 

Global plate models, which yield the relative motion between the North American and 
Pacific plates, have been used as kinematic boundary conditions on the integrated 
deformation across the plate boundary zone in the Western United States [Minster and 
Jordan, 1984, 1987; Bird and Rosenstock, 1984; Weldon and Humphreys, 1986]. If the 
San Andreas fault functioned as a simple boundary that accommodated the full motion 
between two rigid plates, the predicted rate of slip in central California would be ~49 
mm/yr at ~N35°W [DeMets et al., 1987]. Deformation across the Pacific-North American 
boundary is instead thought to be distributed across a broad zone between the continental 
escarpment and the eastern front of the Basin and Range province. As summarized above, 
the rate of slip on the San Andreas fault in central California oriented at N41*W ± 2* is 
constrained from Holocene geological data and ground-based geodetic measurements to be 
34 ± 3 mm/yr. The vector difference between the plate motion and the San Andreas rate is 
-16 mm/yr in the direction N21*W. The integrated rate of deformation to the east of the 
San Andreas fault across the Basin and Range has been estimated from geological 
observations averaged over the Holocene and from Very Long Baseline Interferometry 
(VLBI) to be 9.7 ± 2.1 mm/yr at N56“W ± 10* [Minster and Jordan, 1987]. The vector 
difference derived using the above rates of motion for the San Andreas fault and Basin and 
Range, referred to as the discrepancy vector [Minster and Jordan, 1987], is ~10 mm/yr in 
the direction N14“E or ~5 mm/yr parallel to the San Andreas fault and ~8 mm/yr of 
convergence normal to the fault. On the basis of estimations made by Minster and Jordan 
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[1987], the uncertainities in the discrepancy vector are approximately ±5 mm/yr for the rate 
of slip and ±15* in the direction. Although some minor internal deformation within the 
Sierra Nevada block [Lockwood and Moore , 1979] or across the Great Valley syncline 
may occur, most of the deformation represented by the discrepancy vector is thought to 
occur within the Coast Ranges. 

A compilatation of rates of deformation in the Coast Ranges, separated into estimated 
rates of right-lateral strike-slip motion on specific faults and rates of distributed 
compression, is given in Table 4. In general the geological data are most useful for placing 
upper and lower bounds on rates and for indicating long-term modes of deformation. The 
large range in most geologically determined rate estimates are due to uncertainities in dating 
rock units and in the interpretation of geological reconstructions. To extrapolate 
incremental strain rates determined over a geologically short interval to infer long-term slip 
rates requires that the measured process accumulates strain in a temporally uniform manner. 
Temporal variations in the rate of strain accumulation, however, have been documented to 
occur during and after large earthquakes [Thatcher, 1986], and there is a suggestion that the 
rate of slip on the Paicines and Calaveras faults changed at about the time of the 1979 
Coyote Lake earthquake (see earlier section). An additional complexity involved in using 
geodetic data is that the measured strain rate may be due to the superposition of more than 
one ongoing geologic process. This case is illustrated by geodetic data from the USGS 
San Luis network located near Parkfield, California (Figure 5). Within the area spanned by 
the San Luis network, geological studies suggest two primary modes of neotectonic 
deformation (see discussion of deformation west of the San Andreas fault in the Tectonic 
Setting section). The San Andreas fault in this region ruptures in periodic moderate 
earthquakes on approximately 22-year intervals [Bakun and McEvilly, 1984]; the last 
earthquake occurred in 1966. During the interseismic period, right-lateral shear strain at the 
orientation of the San Andreas fault accumulates. To the west of the San Andreas fault, 
northeast-southwest compression is thought to occur within the Coastal Belt province. 



When subregions of the San Luis network within the Coastal Belt province were examined 
the rate of strain was found to be right-lateral shear strain at the orientation of the San 
Andreas fault [King et al., 1986]. As discussed in Table 4, Harris and Segall [1987] were 
only able to estimate the rate of fault-normal compression by including it as a parameter in 
their inversion. As discussed above, on the segment of the San Andreas fault adjacent to 
the San Benito network, fault slip occurs primarily through creep and we would expect 
little right-lateral shear strain accumulation associated with the San Andreas fault to be 
measured on off-fault geodetic lines. Thus we interpreted the rate of strain determined 
from the San Benito measurements to be associated with the fold structures of the region. 

To explore the implications of the constraints provided by the San Andreas 
discrepancy vector we compare the predicted rates of fault-parallel and fault normal 
deformation to the values given in Table 4. There are several faults on which recent right- 
lateral strike-slip displacements have been observed. At this point we cannot distinguish 
between a model in which 5 mm/yr of predicted fault-parallel motion is distributed over all 
the strike-slip faults in Table 4 and a model in which only the San Gregorio-Hosgri fault 
accommodates the predicted slip rate. The geological and seismicity data (Figure 1) 
suggest that northeast-southwest compression may be localized to two regions, the 30 km- 
wide zone spanned by the San Benito network and within the Coastal Belt province. If 
shortening across the Coast Ranges was divided between these two regions, approximately 
4 mm/yr of shortening is predicted to occur within the Diablo Range. If the 8 mm/yr is 
distributed uniformly across the 170 km- wide zone between the western edge of the Great 
Valley and the continental escarpment, the shortening across the region spanned by the San 
Benito network would be approximately 1 .4 mm/yr. The rate of deformation estimated 
from the San Benito study is most consistent with the first hypothesis; but the uncertainities 
in our calculated values do not allow the latter hypothesis to be ruled out 

Results from ongoing geodetic studies with stations in central California should 
provide additional constraints on the rate and distribution of slip within the Coast Ranges. 
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The rate of slip estimated from VLBI measurements made at Vandenberg, Fort Ord, 
Presidio, and Point Reyes (Figure 5) could potentially be used to constrain the rate of 
deformation across the Coast Ranges. To distinguish between different models will require 
uncertainities in the rate of slip, relative to a fixed reference, of 2-3 mm/yr. Our 
preliminary analysis of VLBI data, based on measurements between October 1982 and 
March 1987, indicate that only the rate of slip at the station Vandenberg meets this 
requirement [Sauber et al., 1987; see also Clark et al., 1987]. Additional measurements of 
all of the VLBI stations in central California have been made within the last year and an 
analysis of these measurements will be the subject of a later paper. Rate of deformation 
estimated from measurements to the Farallon Islands (Figure 5) will be updated following a 
trilateration survey of the network by the USGS in 1988; such information should provide 
a better constraint on the rate of slip on the San Gregorio fault 

Summary 

We have examined triangulation and trilateration data from two geodetic networks that 
lie within the region north of Coalinga, California, between the Great Valley and the San 
Andreas fault The shear strain rate was determined to be 0.15 ± 0.08 prad/yr for the 
Diablo Range region. If we assume that this deformation is due to compression across the 
folds of this area, as is suggested by the orientation of the axis of maximum compressive 
strain (N16*E ±14°), the average shear straining corresponds to a relative shortening rate of 
4.5 ± 2.4 mm/yr. 

We compared the geodetically determined incremental principal strain direction to the 
orientation of maximum principal stress directions estimated from wellbore breakouts and 
the azimuths of P-axes determined from earthquake focal mechanisms 
east of the San Andreas fault in central California. The breakout orientations and the P- 
axes directions within the Diablo Range and near the western edge of the Great Valley were 
found to be similar to the orientation of maximum compressive strain implied by the trend 
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of local fold structures. Our results suggest that the principal stress and strain directions 
within the Diablo Range not support the hypothesis of uniform fault-normal compression 
proposed by Mount and Suppe [1987], 

For a zone within 10 km of the San Andreas fault, trilateration measurements on off- 
fault lines east of the San Andreas fault as well as lines that cross the San Andreas fault 
have been used to estimate the rate of slip along the Calaveras-Paicines fault and to 
document the gradual southward transition in the width of the zone accommodating right- 
lateral fault slip. South of Hollister the inferred rate of slip on the Calaveras-Paicines fault 
was found to be 10-12 mm/yr. The rate of slip on the Paicines fault near Bitter is ~4 
mm/yr. Further south all the slip (at least across the ~20 km wide zone of measurements) 
occurs on the San Andreas fault 
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Figure Captions 


Figure 1. Earthquakes (Ml ^ 4.0) in the Coast Ranges from January 1962 through 
December 1982 [Engdahl and Rinehart, 1987]. The location of three 
earthquakes with well-determined focal mechanisms, the 1982 Idria, 
the 1983 Coalinga, and the 1985 North Kettleman Hills events are given by 
stars. Fault traces are simplified from Jennings [1975]: HF=Hosgri fault, 
RF=Rinconada fault, SAF=San Andreas fault, SGF=San Gregorio fault, 
WTR=Westem Transverse Ranges. SF=City ofSan Fransico, An outline 
of the San Benito triangulation/trilateration network is given for reference. 

Figure 2. Location of stations in the San Benito triangulation/trilateration network 
(solid circles), the Coalinga trilateration network (open circles), and small 
aperture trilateration networks (triangles). Half-filled circles denote stations 
that were part of both the San Benito and Coalinga networks. Quaternary 
faults [Jennings, 1975] are also shown. #C=city of Coalinga, California, 
#H=city of Hollister, California. Stations discussed in the text include 
BIT=Bitter, BRN=Brown, CHL=Chalone, CRS=Cross, HEP=Hepsedam, 
and PAN=Panoche. 

Figure 3. Major structural features of the Diablo Range between Hollister and 

Coalinga [modified from Dibblee, 1979]. The locations of stations in the 
San Benito triangulation/trilateration network are shown as circles for 
reference. #C=city of Coalinga, California, #H=city of Hollister, 
California, CA=Coalinga Anticline, CFZ=Calaveras fault zone, NI=New 
Idria diapir, OFZ=Ortigalita fault zone, PF=Paicines fault, PS=Paicines 
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Syncline, VS=Vallecitos Syncline. Inward pointing double arrows indicate 
a syncline, outward pointing double arrows an anticline. Single arrows 
indicate the direction of plunge of a fold axis. 

Figure 4. Fold structures and measured stress orientation in the region east of the San 
Andreas fault in central California [modified from Mount and Suppe, 1987]. 
Dashed lines indicate a syncline, a dotted line indicates an anticline. Single 
arrows indicate the direction of plunge of a fold axis. Wellbore breakout 
measurements are given by a solid line with inward pointing arrows 
perpendicular to the direction of oi. The direction of Ci inferred from the 
azimuth of P-axis for three earthquakes are labeled 1,2, and 3: 1 = New 
Idria event, 2 = Coalinga event, and 3 = North : Kettleman Hills event. 

The axis of maximum compressive strain ((3) determined from triangulation 
and trilateration data is given along with the uncertainity in the direction. 

Figure 5. Reference figure for the studies given in Table 4. B=Brush and 

M=Mulligan are two stations in the USGS Pajaro Trilateration Network, f 
and g are sites at which geologically estimated rates of fault slip were 
determined for the San Gregorio and Hosgri faults (see f and g in Notes 
for Table 4). The line A-A' is that for which Namson and Davis, [1987] 
constructed their geologic cross section, see i in Notes for Table 4. 
Faults are as given in Figure 1. FIN = Farallon Islands Network, SBN = 
San Benito Network, SLN = San Luis Network. VLBI stations: FTOR = 
Fort Ord, PRES = Presidio, PTRY = Point Reyes, and VNDN = 
Vandenberg. 
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Table l.Strain Rate Parameters for Spatial Subsets of the San Benito 

Network 


• ^ 

Subnet Number of Angles Yl Y2 Y P V 

urad/vr urad/vr urad/vr 


East 

14 

0.18 ±0.11 0.01±0.12 

0.1810.10 

N 1*W±20* 

N46*W±20° 

West 

9 

0.07±0.14 -0.18±0.13 

0.1910.13 

N34*E±22 # 

N1TW122* 

North 

13 

0.04±0.10 -0.12±0.10 

0.1610.11 

N37*E±20* 

N 8°W±20* 

South 

7 

0.22±0.16 -0.13±0.22 

0.2510.20 

N10°E±20 # 

N35*W±20° 

All stations 
except BIT, HEP, 
and PAN 25 

0.13±0.08 -0.0810.08 

0.1510.08* N16’E±14°* 

N29 , W±14 0 ’ 1 


* These strain rate parameters have been scaled by the a posteriori variance factor [Vanicek 
and Krakiwsky, 1986]. 
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Table 2.FauIt-Slip Rates Inferred from Length Changes on Fault-Crossing 

Lines 


Calaveras! Paicines Fault Zone 


Station 1 

Station 2 No. of Obs. Period 

L, 

mm/vr 

Azimuth, 
deg_ _ 

Slip Rate, 
mm/vr 

BRN 

CRS 

7 

72.7-83.9 

-8.7±0.8 

168 

9.9±0.9 

CRS 

CHL 

9 

72.1-83.9 

-18.5±0.4 

181 

24.9±0.5 

BRN 

CHL 

8 

72.1-83.9 

27.4±0.4 

175 

33.8±0.5 


Paicines Fault Adjacent to the Central Creeping Portion of the San Andreas Fault 

• 

Station 1 Station 2 No. of Obs. Period L, Azimuth, 

mm/vr dee 

Slip Rate, 
mm/vr 

BIT 

HEP 

6 

73.1-83.9 

3.7±0.5 

129 

3.8±0.5 

CHL 

BIT 

8 

72.1-83.9 

19.8±0.7 

280 

25.8±0.9 

CHL 

HEP 

7 

72.7-83.9 

27.4±0.6 

113 

30.7±0.6 


L is the rate of line length change determined by least-squares, together with one standard 
deviation. Azimuth is measured clockwise from north. The slip rate is that appropriate 
to the strike-slip fault crossed by the indicated line. 


The time period of observations are given in the decimal fraction of year. 


31 


Table 3.Fault Slip Rates Indicated by Short-Range Trilateration Networks 


Network 

No. of Obs. 

Period 

Fault Slip 
mm/vr 

Pionne 

8 

75.2-87.3 

12±2 

Dry Lake 

4 

79.0-87.3 

27+2 

Tully 

8 

74.9-87.3 

32±1 
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Table 4.Summary of Deformation Rates within the Central Coast Ranges 


Fault 

Orientation 

Geological Slip Rate, 
mm/vr 

Geodetic Slip Rate, 
mm/vr 

Right-lateral strike-slip 



Ortigalita 

N35°W 

0-2 a 


San Andreas 

N41’W 

31-37^ 

32±3 C 

Rinconada 

N35*W 

0-2d 

2±l e 

San Simeon 

N34*W 

6-9^ 


San Gregorio 

N20’W 

7-118 

0±8 h 

Compression in the Coast Ranges 



East of San Andreas 

N49*E 

2.2-5.51 

4.5±2.4 j 




(N16*E) 

West of San Andreas 

N49’E 

4.4-1 li 

6.1±1.7k,l 



• 
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Notes for Table 4 


a Hartetal. [1986]. 
b Siehand Johns [ 1984]. 

c This study; Savage and Burford [1973]; Thatcher [1979b]; Burford and Harsh [1980]; 
Lisowski and Prescott [1981]. 

d Hart et al. [1986]; D. B. Slenunons, personal communication, 1987; E. W. Hart, 
personal communication, 1988. 

e This study. The line between Brush and Mulligan of the Pajaro trilateration network has 
been measured five times between May 1978 and April 1983. If we assume that the 
average line length change is due to right-lateral slip on the King City fault, a northern 
extension of the Rinconada fault, then a slip rate of 2±1 mm/yr is indicated, 
f Clark etal. [1985]. 
g Clark et al. [1985], two sites. 

h On the basis of geodetic measurements made to the Farallon Islands between mid- 1979 
and late \9%5[Prescott and Yu , 1986]. 

i Namson and Davis [1987] obtained 11 km of shortening between the western Great 
Valley and the San Andreas fault. The 22 km of shortening to the west of the San Andreas 
fault was computed by deriving a solution that satisfies the observed structural relief. The 
ranges given for the rate estimates were calculated by assuming that shortening commenced 
between 5 m.y. and 2 m.y. ago. 
j This study. 

k Segall and Harris [1986]; Harris and Segall [1987]. Average rates of line length 
changes from the San Luis trilateration network were used to invert for slip rate at depth on 
the San Andreas fault. In order to fit the trilateration measurements from this network, it 
was necessary to include a component of contraction normal to the trend of the San 
Andreas fault. The inversion results suggest a spatially uniform normal strain of -0.06 
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(istrain/yr. The net shortening rate across the network is 6.1 ± 1.7 mm/yr. This estimated 
compression, however, may be due to a systematic bias in the older trilateration data [J.C. 
Savage, personal communication, 1987 ] . 

12 Two additional geodetic studies using historical triangulation data have been made west 
of the San Andreas fault Burford [1967] analyzed triangulation data measured between 
1930 and 1951 from two networks. One extends from the region where the San Andreas 
and Calaveras faults diverge out to Monterey Bay and one extends from Kettleman Hills 
near Parkfield west to San Luis Obisbo. Outside a zone close to the San Andreas fault, the 
direction of maximum shortening was estimated to be approximately N35*E. As part of a 
general study of the deformation in central California Thatcher [1979b] examined 
triangulation data measured during the time interval 1944-1963 from the Salinias Valley 
network. Most of this network lies within the Salinian block located between the San 
Andreas and Rinconada faults. Examination of subregions suggest that the strains are 
poorly resolved with the orientation of the inferred strain field not correlated with any 
known faults or tectonic trends. Because of the uncertainities in the results, neither of these 
studies were used for rate determinations. 
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